To adapt the use of GH3.TRE-Luc reporter gene cell line for a quantitative high-throughput screening (qHTS) platform, we miniaturized the reporter gene assay to a 1536-well plate format. 1280 chemicals from the Library of Pharmacologically Active Compounds (LOPAC) and the National Toxicology Program (NTP) 1408 compound collection were analyzed to identify potential thyroid hormone receptor (TR) agonists and antagonists. Of the 2688 compounds tested, eight scored as potential TR agonists when the positive hit cut-off was defined at ≥10% efficacy, relative to maximal triiodothyronine (T3) induction, and with only one of those compounds reaching ≥20% efficacy. One common class of compounds positive in the agonist assays were retinoids such as all-trans retinoic acid, which are likely acting via the retinoid-X receptor, the heterodimer partner with the TR. Five potential TR antagonists were identified, including the antiallergy drug tranilast and the anxiolytic drug SB 205384 but also some cytotoxic compounds like 5-fluorouracil. None of the inactive compounds were structurally related to T3, nor had been reported elsewhere to be thyroid hormone disruptors, so false negatives were not detected. None of the low potency (>100µM) TR agonists resembled T3 or T4, thus these may not bind directly in the ligand-binding pocket of the receptor. For TR agonists, in the qHTS, a hit cut-off of ≥20% efficacy at 100 µM may avoid identification of positives with low or no physiological relevance. The miniaturized GH3.TRE-Luc assay offers a promising addition to the in vitro test battery for endocrine disruption, and given the low percentage of compounds testing positive, its high-throughput nature is an important advantage for future toxicological screening.

Endocrine disrupting chemicals (EDCs) interfere with the biosynthesis, metabolism, and functions of steroid hormones, including estrogens and androgens. Aromatase enzyme converts androgen to estrogen. Thus, EDCs against aromatase significantly impact estrogen- and/or androgen-dependent functions, including the development of breast cancer. The current study aimed to develop a biologically relevant cell-based high-throughput screening assay to identify EDCs that act as aromatase inhibitors (AIs), estrogen receptor (ER) agonists, and/or ER antagonists. The AroER tri-screen assay was developed by stable transfection of ER-positive, aromatase-expressing MCF-7 breast cancer cells with an estrogen responsive element (ERE) driven luciferase reporter plasmid. The AroER tri-screen can identify: estrogenic EDCs, which increase luciferase signal without 17β-estradiol (E2); anti-estrogenic EDCs, which inhibit the E2-induced luciferase signal; and AI-like EDCs, which suppress a testosterone-induced luciferase signal. The assay was first optimized in a 96-well plate format and then miniaturized into a 1536-well plate format. The AroER tri-screen was demonstrated to be suitable for high-throughput screening in the 1536-well plate format, with a 6.9-fold signal-to-background ratio, a 5.4% coefficient of variation, and a screening window coefficient (Z-factor) of 0.78. The assay suggested that bisphenol A (BPA) functions mainly as an ER agonist. Results from screening the 446 drugs in the National Institutes of Health Clinical Collection revealed 106 compounds that modulated ER and/or aromatase activities. Among these, two AIs (bifonazole and oxiconazole) and one ER agonist (paroxetine) were confirmed through alternative aromatase and ER activity assays. These findings indicate that AroER tri-screen is a useful high-throughput screening system for identifying ER ligands and aromatase-inhibiting chemicals.

Chordoma is a rare primary bone malignancy that arises in the skull base, spine and sacrum and originates from remnants of the notochord. These tumors are typically resistant to conventional chemotherapy, and to date there are no FDA-approved agents to treat chordoma. The lack of in vivo models of chordoma has impeded the development of new therapies for this tumor. Primary tumor from a sacral chordoma was xenografted into NOD/SCID/IL-2R γ-null mice. The xenograft is serially transplantable and was characterized by both gene expression analysis and whole genome SNP genotyping. The NIH Chemical Genomics Center performed high-throughput screening of 2,816 compounds using two established chordoma cell lines, U-CH1 and U-CH2B. The screen yielded several compounds that showed activity and two, sunitinib and bortezomib, were tested in the xenograft. Both agents slowed the growth of the xenograft tumor. Sensitivity to an inhibitor of IκB, as well as inhibition of an NF-κB gene expression signature demonstrated the importance of NF-κB signaling for chordoma growth. This serially transplantable chordoma xenograft is thus a practical model to study chordomas and perform in vivo preclinical drug testing.

Chronic exposure to drinking water disinfection byproducts has been linked to adverse health risks. The monohaloacetic acids (monoHAAs) are generated as byproducts during the disinfection of drinking water and are cytotoxic, genotoxic, mutagenic, and teratogenic. Iodoacetic acid toxicity was mitigated by antioxidants, suggesting the involvement of oxidative stress. Other monoHAAs may share a similar mode of action. Each monoHAA generated a significant concentration-response increase in the expression of a β-lactamase reporter under the control of the antioxidant response element (ARE). The monoHAAs generated oxidative stress with a rank order of iodoacetic acid (IAA) > bromoacetic acid (BAA) ≫ chloroacetic acid (CAA); this rank order was observed with other toxicological end points. Toxicogenomic analysis was conducted with a nontransformed human intestinal epithelial cell line (FHs 74 Int). Exposure to the monoHAAs altered the transcription levels of multiple oxidative stress responsive genes, indicating that each exposure generated oxidative stress. The transcriptome profiles showed an increase in thioredoxin reductase 1 (TXNRD1) and sulfiredoxin (SRXN1), suggesting peroxiredoxin proteins had been oxidized during monoHAA exposures. Three possible sources of reactive oxygen species were identified, the hypohalous acid generating peroxidase enzymes lactoperoxidase (LPO) and myeloperoxidase (MPO), nicotinamide adenine dinucleotide phosphate (NADPH)-dependent oxidase 5 (NOX5), and PTGS2 (COX-2) mediated arachidonic acid metabolism. Each monoHAA exposure caused an increase in COX-2 mRNA levels. These data provide a functional association between monoHAA exposure and adverse health outcomes such as oxidative stress, inflammation, and cancer.

Drug-induced phospholipidosis is characterized by the accumulation of intracellular phospholipids in cells exposed to cationic amphiphilic drugs. The appearance of unicentric or multicentric multilamellar bodies viewed under an electron microscope (EM) is the morphological hallmark of phospholipidosis. Although the EM method is the gold standard for detecting cellular phospholipidosis, this method has its drawbacks, including low throughput, high cost, and unsuitability for screening a large chemical library. In this study, a cell-based phospholipidosis assay has been developed using the LipidTOX Red reagent in HepG2 cells and miniaturized into a 1536-well plate format. To validate this assay for high-throughput screening (HTS), the LOPAC library of 1280 compounds was screened using a quantitative HTS platform. A group of known phospholipidosis inducers, such as amiodarone, propranolol, chlorpromazine, desipramine, promazine, clomipramine, and amitriptyline, was identified by the screen, consistent with previous reports. Several novel phospholipidosis inducers, including NAN-190, ebastine, GR127935, and cis-(Z)-flupentixol, were identified in this study and confirmed using the EM method. These results demonstrate that this assay can be used to evaluate and profile large numbers of chemicals for drug-induced phospholipidosis.

A goal of the Tox21 program is to transit toxicity testing from traditional in vivo models to in vitro assays that assess how chemicals affect cellular responses and toxicity pathways. A critical contribution of the NIH Chemical Genomics center (NCGC) to the Tox21 program is the implementation of a quantitative high throughput screening (qHTS) approach, using cell- and biochemical-based assays to generate toxicological profiles for thousands of environmental compounds. Here, we evaluated the effect of chemical compounds on mitochondrial membrane potential in HepG2 cells by screening a library of 1,408 compounds provided by the National Toxicology Program (NTP) in a qHTS platform. Compounds were screened over 14 concentrations, and results showed that 91 and 88 compounds disrupted mitochondrial membrane potential after treatment for 1 or 5 h, respectively. Seventy-six compounds active at both time points were clustered by structural similarity, producing 11 clusters and 23 singletons. Thirty-eight compounds covering most of the active chemical space were more extensively evaluated. Thirty-six of the 38 compounds were confirmed to disrupt mitochondrial membrane potential using a fluorescence plate reader, and 35 were confirmed using a high content imaging approach. Among the 38 compounds, 4 and 6 induced LDH release, a measure of cytotoxicity, at 1 or 5 h, respectively. Compounds were further assessed for mechanism of action (MOA) by measuring changes in oxygen consumption rate, which enabled the identification of 20 compounds as uncouplers. This comprehensive approach allows for the evaluation of thousands of environmental chemicals for mitochondrial toxicity and identification of possible MOAs.

Both pharmacophore models of the human ether-à-go-go-related gene (hERG) channel blockers and phospholipidosis (PLD) inducers contain a hydrophobic moiety and a hydrophilic motif/positively charged center, so it is interesting to investigate the overlap between the ligand chemical spaces of both targets. We have assayed over 4000 non-redundant drug-like compounds for both their hERG inhibitory activity and PLD inducing potential in a quantitative high throughput screening (qHTS) format. Seventy-seven percent of PLD inducing compounds identified from the screening were also found to be hERG channel blockers, and 96.9% of the dually active compounds were positively charged. Among the 48 compounds that induced PLD without inhibiting hERG channel, 24 compounds (50.0%) carried steroidal structures. According to our results, hERG channel blockers and PLD inducers share a large chemical space. In addition, a positively charged hERG channel blocker will most likely induce PLD, while a steroid PLD inducer is less likely a hERG channel blocker.

Since its establishment in 2008, the US Tox21 inter-agency collaboration has made great progress in developing and evaluating cellular models for the evaluation of environmental chemicals as a proof of principle. Currently, the program has entered its production phase (Tox21 Phase II) focusing initially on the areas of modulation of nuclear receptors and stress response pathways. During Tox21 Phase II, the set of chemicals to be tested has been expanded to nearly 10,000 (10K) compounds and a fully automated screening platform has been implemented. The Tox21 robotic system combined with informatics efforts is capable of screening and profiling the collection of 10K environmental chemicals in triplicate in a week. In this article, we describe the Tox21 screening process, compound library preparation, data processing, and robotic system validation.

Bisphenol A (BPA), is a well-known endocrine disruptor compound (EDC) that affects the normal development and function of the female and male reproductive system, however the mechanisms of action remain unclear. To investigate the molecular mechanisms of how BPA may affect ten different nuclear receptors, stable cell lines containing individual nuclear receptor ligand binding domain (LBD)-linked to the β-Gal reporter were examined by a quantitative high throughput screening (qHTS) format in the Tox21 Screening Program of the NIH. The results showed that two receptors, estrogen receptor alpha (ERα) and androgen receptor (AR), are affected by BPA in opposite direction. To confirm the observed effects of BPA on ERα and AR, we performed transient transfection experiments with full-length receptors and their corresponding response elements linked to luciferase reporters. We also included in this study two BPA analogs, bisphenol AF (BPAF) and bisphenol S (BPS). As seen in African green monkey kidney CV1 cells, the present study confirmed that BPA and BPAF act as ERα agonists (half maximal effective concentration EC50 of 10-100 nM) and as AR antagonists (half maximal inhibitory concentration IC50 of 1-2 μM). Both BPA and BPAF antagonized AR function via competitive inhibition of the action of synthetic androgen R1881. BPS with lower estrogenic activity (EC50 of 2.2 μM), did not compete with R1881 for AR binding, when tested at 30 μM. Finally, the effects of BPA were also evaluated in a nuclear translocation assays using EGPF-tagged receptors. Similar to 17β-estradiol (E2) which was used as control, BPA was able to enhance ERα nuclear foci formation but at a 100-fold higher concentration. Although BPA was able to bind AR, the nuclear translocation was reduced. Furthermore, BPA was unable to induce functional foci in the nuclei and is consistent with the transient transfection study that BPA is unable to activate AR.

The human cytochrome P450 (CYP) enzyme family is involved in the biotransformation of many xenobiotics. As part of the U.S. Tox21 Phase I effort, we profiled the CYP activity of approximately three thousand compounds, primarily those of environmental concern, against human CYP1A2, CYP2C19, CYP2C9, CYP2D6, and CYP3A4 isoforms in a quantitative high throughput screening (qHTS) format. In order to evaluate the extent to which computational models built from a drug-like library screened in these five CYP assays under the same conditions can accurately predict the outcome of an environmental compound library, five support vector machines (SVM) models built from over 17,000 drug-like compounds were challenged to predict the CYP activities of the Tox21 compound collection. Although a large fraction of the test compounds fall outside of the applicability domain (AD) of the models, as measured by k-nearest neighbor (k-NN) similarities, the predictions were largely accurate for CYP1A2, CYP2C9, and CYP3A4 ioszymes with area under the receiver operator characteristic curves (AUC-ROC) ranging between 0.82 and 0.84. The lower predictive power of the CYP2C19 model (AUC-ROC = 0.76) is caused by experimental errors and that of the CYP2D6 model (AUC-ROC = 0.76) can be rescued by rebalancing the training data. Our results demonstrate that decomposing molecules into atom types enhanced the coverage of the AD and that computational models built from drug-like molecules can be used to predict the ability of non-drug like compounds to interact with these CYPs.